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Abstract 


In this paper we have proposed an almost unbiased estimator using known value of some 
population parameter(s). Various existing estimators are shown particular members of the 
proposed estimator. Under simple random sampling without replacement (SRSWOR) scheme the 
expressions for bias and mean square error (MSE) are derived. The study is extended to the two 
phase sampling. Empirical study is carried out to demonstrate the superiority of the proposed 
estimator. 
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1. Introduction 


Consider a finite population U = Uj, Us, ...Uy of N units. Let y and x stand for the 
variable under study and auxiliary variable respectively. Let (yj, xj), I=1, 2,.., n 
denote the values of the units included in a sample s, of size n drawn by simple 
random sampling without replacement (GSRSWOR). The auxiliary information has 
been used in improving the precision of the estimate of a parameter (See 
Cochran (1977), Sukhatme et. al. (1984) and the references cited there in). Out 
of many methods, ratio and product methods of estimation are good illustrations 
in this context. 


In order to have a survey estimate of the population mean Y of the study 
character y, assuming the knowledge of the population mean X of the auxiliary 
character x, the well-known ratio estimator is 


(1.1) 


Bahl and Tuteja (1991) suggested an exponential ratio type estimator as — 


tan = V CX X =x 
re y Xa 
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several authors have used prior value of certain population parameter(s) to find 
more precise estimates. Sisodiya and Dwivedi (1981), Sen (1978) and 
Upadhyaya and Singh (1984) used the known coefficient of variation (CV) of the 
auxiliary character for estimating population mean of a study character in ratio 
method of estimation. The use of prior value of coefficient of kurtosis in 
estimating the population variance of study character y was first made by Singh 
et. al. (1973). Later used by Singh and Kakaran (1993) in the estimation of 
population mean of study character. Singh and Tailor (2003) proposed a modified 
ratio estimator by using the known value of correlation coefficient. Kadilar and 
Cingi (2006), Khosnevisan et. al. (2007), Singh et. al. (2007) Singh and Kumar 
(2009) and Singh et. al. (2009) have suggested modified ratio estimators by 
using different pairs of known value of population parameter(s). 


In this paper under SRSWOR, we have proposed almost unbiased estimator for 
estimating Y. 
2. Almost unbiased ratio type estimator 


Suppose 
| _ | _ (aX+b — _ ((aX+b)-(ax+b 
fo =, ts = 7( ), Use = YEXP (expec 





x+b 


such that ty t,,,t,..€ Ww, where w, denotes the set of all possible ratio type 
estimators for estimating the population mean Y. By definition the set w, is a 
linear variety, if 

tor — Wo + 0,T.. + WTrse ew (2.1) 
for 3250, =1, w,eR (2.2) 
where w, (i=0, 1, 2) denotes the statistical constants and R denotes the set of 
real numbers. 


To obtain the bias and MSE of t.,,, we write 
y=Y(l+e,), x=X(l+e,), 











such that 

E (eo) =E (e;)=0. 
E(e9) = £\Cy, E(e7) =1,C2, E(ege,) = f)pCyCx 

1 1 1 J —_\p 1 w —\y 

here f, =(—-—), ee -Y), St= -xXxy), 
where fy =(—-—) way woe ) 

S S C S 
C. =, C= K=p| — |, p=7——V, 
yy eX fee PBS.) 
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Expressing t,,, in terms of e’s, we have 


ty =¥(1 +e) |W + ,(1 + 0e,)"1 + wexp {— + (1+ 6e,)7*}] (2.3) 
aX 
where §@ = 
aX+b 


Expanding the right hand side of (2.3) and retaining terms up to second order of 
e's, we have 


t,, =Y E +e, —w6e, + 67 (wo, + =) e. — Bwere, | (2.4) 
where w= (w, +=). 

Taking expectations of both side of (2.4) and then subtracting Y from both side, 
we get the bias of the estimator t,,, up to the first order of approximation as 








B(t,,)=£,Y je*c2 (wo, ~<a ) = Gwe, C,| (2.6) 
B(t,,) = f, ¥[6*c2 — Opc,c, | (2.7) 
26° 5 @pc..Cc., 
B(t ee f,Y = =] (2.8) 
From (2.4), we have 
(t,. —Y¥) = Y[e, — Owe, | (2.9) 


Squaring both sides of (2.9) and then taking expectations, we get MSE of the 
estimator t,,, up to the first order of approximation, as 


MSE(t,,) = f,¥ [CZ + 6700°C2 — 28wpC,C,| (2.10) 
This is minimum when 

w= K(=pZ). (2.11) 
Putting this value of (= k) in (2.10), we get the minimum MSE of t,, as 

min. MSE(t,,) = f, ¥7C2(1— p?) (2.12) 


which is same as that of traditional linear regression estimator 
from (2.5) and (2.11), we have 


@, = (2.13) 


From (2.2) and (2.13), we have only two equations in three unknowns. It is not 
possible to find the unique values for w,;..,i = 0,1,2. In order to get unique values 


for «..., we shall impose the linear restriction 
w)B(¥) + w,B(t,,.) + w,B(t,,.) = 0 (2.14) 
a (2. ot (2.11) and (2.14) can be written as in the matrix form as 
1 1 
lo: a is if (2.15) 
0 B(t,.) Bie.) 


Pak.j.stat.oper.res. Vol.IV No.2 2008 pp63-76 65 


Rajesh Singh, Mukesh Kumar, Florentin Smarandache 


Using (2.15) ,we get unique values of w,., (i=0,1,2) as 

Ag 

Ay 
moe 
where 

A,> a(t) —= B(t,.) 
_ poe 1358 a he 

A.o= B(t,.. {1 ~ k} + 2 B(t,..) {k | | (2.1 7) 


A,= k.B(t,,..) 
A.,= —k B(t,,) 


Use of these «,:, (i=0,1,2) remove the bias up to terms of order (n’') at (2.1). 


3. Product —type estimators 


ciated | _ ( (a¥+b)—(aX+b) 
PHPpOse hs =9(= tose = aed 





such that to, the, tose ©Q, where Q denotes the set of all possible product —type 
estimators for estimating the population mean Y. By definition, the set Q Is linear 
variety if 


t, = doy 7 ditys + Gatyce EQ (3.1) 
for 32.,4;=1, g,€R (3.2) 


where q,(i=0,1,2) denotes the statistical constants. 


Expressing t, in terms of e’s, we have 


ty =F(1+ep) [ooo + w;(4 + Bey) + w,exp (4 (4 + Be,)~}| (3.3) 
where 6 = =, 
aX+b 


Expanding the right hand side of (3.3) and retaining terms up to second power of 
e's, we have 


t, ¥¥[1+e, +Oqe, Ler + qBere, | (3.4) 
where q= q, + 2 (3.5) 
Taking expectations of both sides of (3.4) and then subtracting Y from both 


sides, we get the bias of the estimator t,, up to the first order of approximation 
as 


B(t,) = f,¥ [-= a°C2 + q@pC,C ul (3.6) 
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Bias expression for the estimators t,. and t,.. iS given by 
B(t,.) =f, Y [@pC,c, | (3.7) 

= 1»? Bp Ow 
B(to..) = f,¥ |-ze c2 = (3.8) 
From (3.4), we have 
(t, —Y) = Ye) + Oge,] (3.9) 


Squaring both the sides of (3.9) and then taking expectations, we get MSE of the 
estimator t,, up to the first order of approximation, as 


MSE(t, ) = £,¥? [c2 + 62q2c2 + 26qpC, C, ] (3.10) 
which is minimum for 

q= -k=-p2 (3.11) 
Putting this value of q(=-k) in (3.10),we get the minimum MSE of t, as 

min. MSE(t,) = f,¥°C2(1 — p?) (3.12) 
which is same as that of traditional linear regression estimator. 


From (3.5) and (3.11), we have 
i ieee (3.13) 


From (3.2) and (3.13), we have only two equations in three unknowns. It is not 
possible to find the unique values for qs, i=0,1,2. In order to get unique values of 
qi s, we shall impose the linear restriction 


q0B(¥) 7 qB(t,.) ur q,B(t,..) =0 (3.1 4) 


- . (3.2),(3.13) and (3.14) can be written in the matrix form as 


1 
0 ‘i lk ‘HT (3.15) 
0 B(ty.) B(tyse) 


Solving (3.15), we get the unique values of qi's (i=0,1,2) as- 


Apo 
do = re 
q3 = = (3.16) 
_ Ap2 
3 An 
where 
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= B(t,..) —= B(t,.) 
Ano= mare +k}+ B(tp.){—k— 5} 
Ani= —k Blt) 

A,2=k B(t,.) 


(3.17) 





Use of these q’’s (i=0,1,2) remove the bias up to terms of order o(n’') at (3.1). 


In Appendix A we have listed some of the important known estimators of the 
population mean, which can be obtained by suitable choice of constants 
w,,i=0,1,2, q,,i=0,1,2 andaandb. 


4. Proposed estimators in two phase sampling 


When X is unknown, it is sometimes estimated from a preliminary large sample of 


size n’ on which only the characteristic x is measured (for details see Singh et. al. 
(2007)). Then a second phase sample of size n (n <n) is drawn on which both y 


and x characteristics are measured. Let x = — yr -, x, denote the sample mean of 
x based on first phase sample of size n’, ¥ = => -,y; and x= ~ Yr Xp be the 
Sample means of y and x respectively based on second phase of size n. 

In two phase sampling the estimator t,,. will take the following form 

twa — Wo yt Wid trsd z 5 Mog Crsed eWay (4.1 ) 
for 24 = 1, w,4€R (4.2) 


(aX +b) (ax+b 
where ty; = y(=*) and trse = yex exp [esa 


To obtain the bias and MSE of t,,,,, we write 
y= Y(l+e,), x=X(l+e,), Z =K(it+e,) 
such that 

E (0) = E (e1) = E(e1) =0. 


E(eg) =f Cy, E(ep)=f;Cx,  E(e't)=f2Cy  E(ege,) =f pC Cx 
E(ege') ) = fopC Cx E(eye, ) = BC 

1 1 1 1 
where ty ae ” ~ Cy? 


Following the procedure mentioned in section 2 and 3, we get bias and MSE of 
twa as 


B(tua) = y [0°C2f, (wg + 8) — Op Cy C.F (wy + 24)| (4.3) 
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MSE(t..4) = Y"[£,C2 + f,C2w4— 2f pC,C, wa] (4.4) 

1 1 
where fz = -——=f, — fy. 

n n 
MSE (t,,4 ) iS minimum, when 

wd 

Wig t=>— == k. (4.5) 
Putting this value of w, in (4.4), we get the minimum MSE of t,,; as 
min. MSE(t.,, ) = ¥°C2[f, —f,p7] (4.6) 
This is same as that of traditional two phase linear regression estimator. 
The bias expression for the estimators t,., and t,.3. is respectively given by 
B(tyeq) = ¥[87C2f,w14 — OpC,C,f,0044| (4.7) 
B(treae) = ¥ G c2f, 2 — 800, C.f, 8] (4.8) 


From (4.2) and (4.5), we have only two equations in three unknowns. It is not 
possible to find the unique values for Wigs I=0,1,2. 


In order to get unique values of w;,, we shall impose linear restriction 








Wog B(¥) - @,4B(t,,4) + 54 B(teea) =0 (4.9) 
vanes (4.2), (4.5) and (4.9) can be written in matrix form, as 
1 
0 a | ose} ; (4.10) 
0 a(t...) B(t,.e4) 
Solving (4.10), we get the unique values of w,.s, (i = 0,1,2) as 
Noa 
Onq = — 
od Ara 
_ Asa | 
Oa = 5 (4.11) 
O24 = = 
ra 
where 
Ava= B (treoa) ~ ; B (trea) 
| | 1 7 1 
Avoa= B(treea) {1 7 k} + 2 B(t,.a) {k- =} 3 (4.12) 


Aua= k. B(t,.ea) 

Ava —k .B(t,.3) 
Use of these wia’s (i=0,1,2) will remove the bias up to terms of order O(n") at 
(4.1). 
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The estimator t, written in (3.1), in two phase sampling, will take following form 
tog = Goa¥ + Giatpsad + d2atpsead € Qa (4.13) 
For di-odia =1, aa eR (4.14) 
where q;,(i=0,1,2) denotes the statistical constants . 
The estimators t,.3 and t,.ge are 

_ ( ax+b 
tosa = 7(=—) and 


t {== =| 
psed — Yexp (ax+b)+(ax'+ 


Following the procedure of section 4, we get the unique values of digs (i=0,1,2) 











as 

—  _ Apod 

od - 

qia= yy (4.15) 
_ Ap2a 

2d ray 

where 


A,a= B(teed) 7 : B( toa) 
Anoa= B(traea ){1 + k} + E B(t,.a) {—k 7 +} | (4.16) 
Apia —k B(treea ) 
Anza =k.B (trea ) 


where 
B(t,.4) - ¥[@q,4f3eC,,C.,| (4.17) 
B(tosea) = ¥[0 M2 f,9C,C, —< arg fC?| (4.18) 


The minimum MSE of t,, is given by 
MSE(t,4) = Y°C>[f, — f,p°) . 


5. Empirical study 


For empirical study we use the data sets earlier used by Kadilar and Cingi (2006) 
(population 1) and Khosnevisan et. al. (2007) (population 2) to verify the 
theoretical results. 
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Data statistics 
TTS 77 6 [pw 
Population? | 20 | 8 | 1965 | 188 | 0.956 | 0004 | -o92 | 3.06 





Table 5.1: Values of w;’'s and gis 


pws Population 1 Population 2 






8.590718 7.892148 
q (21.417) ( 2.919085) 
©, 11.86615 5.234461 
‘ch (16.14158) ( 3.576773) 
Ww, -19.4569 -12.1266 
(q2) ( -36.5586) ( -5.49586) 
(PRE) 


The percent relative efficiencies (PRE) of various estimators of Y are computed 
and presented in Table 5.2 below. 


Table 5.2: PRE of different estimators of Y 


| Estimator | PRE (Pop!) | Estimator | PRE (Pop Il) _ 
sto | 100 | | 100 
pts | 212810 | gs | 465.501 
| ts | 143.923 | gs | 285.920 
| to | 143.988 | gs | 338.487 
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In order to see the performance of the suggested estimators in two phase 
Sampling we use the data set of Murthy (1967) (Population III) and Steel and 
Torrie (1960) (Population IV). 


Population | c, ] G |p | N_ 


|Population1 —_| 0.3542 | 0.9484 | 0.9150 | 80 | 30 | 10 | 
Population2 0.4803 | 0.7493 | -0.4996 


Table 5.3 : The values of w,,‘s and q;,.'s 





Population Ill Population IV 
Woa -0.241523 3.011435 
Toa) (1.808833) (1.089979) 
Wig -0.558071 1.370950 
(Q1a) (0.125381) 0.730464) 


1.799595 -3.382385 
(-0.934214) (-0.820443 ) 





The percent relative efficiencies of various estimators of Y in two phase sampling 
are computed and presented in Table 5.4 below. 


Table 5.4 : PRE of different estimators of Y in two phase sampling 


| Estimator | PRE (Population |) | PRE (Population Il) 
36.642 24.562 
4.849 59.770 
200.420 48.365 

t 


23.628 115.142 
276.156 63.452 


123.762 
123.762 





6. Conclusion 


From theoretical discussion and empirical study we conclude that the proposed 
estimators under optimum conditions perform better than other estimators 
considered in the article. The relative efficiencies of various estimators are listed 
In Table 5.2 and 5.4. 
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Appendix A 


Table A.1: Some members of the proposed family of estimators - 


Ratio Estimator 


(corresponding to 
@ ; i=0,1 52) 


The mean per unit estimator 


aA 
t. =V- 
i~“Vz 


The usual ratio estimator 


x+C, 


Sisodia and Dwivedi (1981) 
estimator 


t; 7 X +B) 


Singh et. al. (2004) estimator 
_— XB2(x)+Cx 
XB (x) + Cy, 
Upadhyaya and Singh (1999) 
estimator 
_— XC, +B2(x) 
XC, +B (x) 


Upadhyaya and Singh (1999) 
estimator 


 _ _|X+p 
t, =y Zio 


Singh and Tailor (2003) 
estimator 


t7 =yex X =x 
aa Al ee 


Bahl and Tuteja (1991) 


estimator 
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Product Estimator 
(corresponding to q;.i=0,1,2) 


Go ¥ 


The mean per unit estimator 
_ 

qi~ VR 

The usual product estimator 

2 aes 

X44, 


Pandey and Dubey (1988) 
estimator 


q2 


_X+B5(x) 
q3 => 
X +B (x) 

Singh et. al. (2004) estimator 
_s XB (x) TF Ce 

XBo(x)+C, 
Upadhyaya and Singh (1999) 
estimator 
qs = _ xC, + B(x) 

XC, F B(x) 


Upadhyaya and Singh (1999) 
estimator 


“ofr 


Singh and Tailor (2003) 
estimator 


= yex X=X 
q7 y MX 


W4 


Bahl and Tuteja (1991) 


estimator 
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Ratio Estimator 


(corresponding to 
i=0,1,2) 


@ 5 Jl 


A 
9 wPEX 45 42B,(0) 


Singh et. al. (2007) estimator 


to = vex _ X-x 
o—# pees eren 


Singh et.al. (2007) estimator 


tin = yexp| = 
10-—Y Wes 


Singh et. al.(2007) estimator 


Sex OK 
‘1 BOO K4AN)42C 
Singh et. al. (2007) estimator 


1 —jexp— 2 ™ 
C(K49)42800) 


Singh et. al. (2007) estimator 


_ C(X—x) 
13=Yyexp ——= 
C(X+x)+20 


Singh et. al. (2007) estimator 


Singh et. al. (2007) estimator 


BOX -X) 
B@X+x)+20 
Singh et. al. (2007) estimator 


Le=y ex} A 
[AX +X)4+-28, OO) 


Singh et. al. (2007) estimator 


t,5=yex 


Product Estimator 


(corresponding to q,,i=0,1,2) 


: sexy 
X+X+26, (x) 


Singh et. al. (2007) estimator 


ae ont A 2C. 

Singh et.al. (2007) estimator 
exp| —— as 

dio = YEXP ene 2p 


Singh et. al.(2007) estimator 


BOK-X) 
BOVK+9)+2G 
Singh et. al. (2007) estimator 


 [ G@X _ 
a5 ae C (XX) +2300 


singh et. al. (2007) estimator 


_ C.(x—X) 
413 ~ PC (X+5) +20 


qd, s=yex 


singh et. al. (2007) estimator 


p(x—X) | 


Cx 
“y .=yerp| Pe 


singh et. al. (2007) estimator 


8,2) @—X) | 
BO)(X+X)+20 
Singh et. al. (2007) estimator 


_ pe—®) 
he“ YON TX 428,00) 


Singh et. al. (2007) estimator 


Gis =Yex 
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In addition to above estimators a large number of estimators can also be 
generated from the proposed estimators just by putting different values of 
constants w-,i=0,1,2, q;,i=0,1,2 andaandb. 
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